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I. INTRODUCTION 

Previous studies of the thermal conductivity of lead 
crystals in the superconducting state, which were de- 
formed plastically by low temperature stretching of the 
initially perfect samples, and observation of the recovery 
processes on annealing of the samples at room tempera- 
tures, had demonstrated strong anomalies in the thermal 
conductivity of the deformed Pb crystals below 4 K.— The 
same effects were also seen in weakly bent Bi crystals? 2 - 
Furthermore, experiments on the thermal conductivity 
of hep 4 Hc crystals grown from high pure 4 He in a long 
capillary had also revealed strong anomalies in thermal 
conductivity of samples that were weakly deformed by 
bending them at temperatures near and above 0.4 K££ 

Several attempts for a theoretical explanation of these 
results have been made, but none have unfortunately 
been completely succesfull.— In this primer paper how- 
ever, we introduce a new model for explaning the ob- 
served anomalies in the thermal conductivity of the 
weakly deformed crystals from high pure matter. This 
model is based on phonon scattering on mobile kinks on 
the newly induced dislocation lines. Previously, a similar 
model, based on scattering of electrons by mobile kinks, 
has been introduced for the explanation of the anomaly 
in the electronic contribution to the thermal conductiv- 
ity of plastically deformed copper crystals^ In systems 
where the phonon thermal conductivity is the main con- 
tribution to the transfer of heat flux, such as quantum 
crystals, metal crystals in superconducting state and non- 
metals, the scattering of thermal phonons by the mobile 
kinks on dislocation lines induced under weak deforma- 
tion of initially perfect samples at reduced temperatures 
seems to be the natural explanation of the experimen- 
tally observed effects. This paper will only introduce 
this process and show the main results of detailed calcu- 
lations of the thermal conductivity in different directions 
relative to the glide plane of the dislocations. We have 
found that in the crystals where scattering of phonons on 



kinks is the dominant scattering process our theoretical 
results can qualitatively reproduce the experimental fea- 
tures. The detailed calculations referred to in this primer 
note and the quantitative fit of the experimental results 
can be found in a paper which is soon to appear J- 

II. KINEMATICS 

For a description of the kinematics of phonon-kink 
scattering we use a similar procedure from Ref. H|. We 
consider a crystal which contains dislocations due to an 
external influence on the crystal. The dislocations lie in 
the xz plane and the direction parallel to the dislocations 
is the z direction. 

Around a dislocation the displacement Uj can be de- 
composed in two components 

u j =u' j +uf. (2.1) 

The " static" displacement itj depends on the presence of 
the kinks and can be written as 

u i = E &( r -L : ^U^V K[z - Mt)) + u' j0 , (2.2) 

K 

where £o( K ) is the Fourier transform of the dislocation's 
line displacement due to the kink, fj(j± : n) is a pro- 
portionality constant and is displacement around the 
straight dislocation without kink. The abbreviation rj_ 
indicates {x,y). The "dynamical" displacement uj has 
its origin in the phonons and can be expressed as a su- 
perposition of plane waves, 

^ = ^g(k, S ) ej (k, S K k - r , (2.3) 

k,s 

where s indicates the polarization of the lattice vibrations 
and e the polarization vector. Treating the kink in a 
harmonic trap (potential well) with angular frequency 
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n and "writing wo(k, s) for the angular frequency of the 
phonons results in the total Lagrangian 



L = L q + L Z(1 + Li n t + const. 
£, = ^£{<KM)«*(M 



k..s 



-Wo(k,s)g(k,s)<f (k,s)} , 
M 



(2.4a) 



(2.4b) 



L Z0 = Y {m - n 2 (z a (t) - z° ) 2 } , (2.4c) 
A„t = -ipV kMt)e-^ Zo{t) Uk z ) 

xF j (k){'(k )> )eJ(k )J ). (2.4d) 

In the equations above, p is the density of the crystal, 
V — L 3 its total volume, M is the kink mass^, zq in- 
dicates the position of the kink, Zq is its rest position 
and Fj(k) is the Fourier transform of fj(r± : k), being 
defined as 

F,(k) = -LJ t < k '• / ; ;r. : k z )d 2 r x , (2.5) 

where k^ = (k x ,k y ). 

From the interaction term L- mt one can determine 
the phonon-kink scattering amplitude per unit time 
A(k, s;k', s'). Due to phonon-kink scattering, phonons 
are no longer described by the Bose-Einstein distribu- 
tion iV (wo(k, s)). In the presence of a small temperature 
gradient VT, the linear correction to the Bose-Einstein 
distribution 5Nk s is given by 

- ^fi AT°Mk, .))(! + iv° M k, ,)) VT ■ 

= E / |^3^(k ) *;k , y)[«JVk i -^ kV ], (2.6) 



with 
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iv°Mk, s ))(i + iv°K(k, s )))- 



(2.7) 



and 



7>(k, S ;k',s') = N ph L 2 \A(k,s;k',s')\ 2 
x lv(w (k, s) - w (k',s');fe,<7z) 
x7VO (wo(k ' ;S ' ))( i +iV K(k, s ))), (2.8) 

with N p h the number of phonons in the crystal and 

K(u; q x , q z ) = — J dzdz' dt exp[iq z (z' — z) + iuit] 
< exp[-iq x £(z, 0)] exp[ig x £(z', t)] > . (2.9) 



With Eq. (|2.6|) a full kinematical treatment of the 
phonon-kink scattering is possible. 



III. HEAT FLOW 

With the full kinematics of the phonon-kink scattering 
at our disposal we arc able to study the effect of phono- 
kink scattering on the heat flow through the crystal. The 
heat flux Q is given by 



Q 



d 3 k aw (k,s) 

fiw (k, s) - diVks 



(27T) 



9k 



-XVT, 



(3-1) 

where x is the matrix of the thermal conductivity. For 
simplicity, we will assume here that this matrix only has 
two distinct diagonal elements and no off-diagonal ele- 
ments 



(3.2) 



This implies that there two distinct heat flows. One along 
the dislocation, 




Q\\=- 

and one perpendicular to, 



-X||(VT) 2 



Qx = -X±(VT) X , 



(3.3) 



(3.4) 



with (VT)_l = ((VT) x ,{VT) y ,0). 

Combined Eqs. (|2.6I) and (|3.ip allow for a full calcu- 
lation^ of x\\ an d X±- This full calculation shows that 
there arc four different temperature regimes for the ther- 
mal conductivity. These four intervals arc 



regime 1: T <C T u , 
regime 2: r u « T « T Q , 
regime 3: T < T < T*, 
regime 4: T > T*, 



Here. 



Tn = - — , 

k B 



2Ain 2 £ 2 



k. 



(3.5a) 
(3.5b) 
(3.5c) 
(3.5d) 
(3.5e) 



(3.6a) 
(3.6b) 
(3.6c) 



where I is the typical size of the kink and ujo(1/£) is 
the angular frequency for a phonon with a wavelength 
equal to the size of the kink. The three temperatures are 
ordered as follows 



T u < To < T* 



(3.7) 



In the calculations we also took into account that in 
real experiments, one does not measure the thermal con- 
ductivity in one particular direction, but rather an aver- 
age over different direction as one has no perfect control 
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of the orientation of the kinks. As the scattering in dif- 
ferent directions is a consecutive process, the scattering 
rates for the different processes add. This means that the 
measured thermal conductivity x is found from 



and at the highest temperatures as 



r 1 



PXI 1 + (1 - P)Xn 



(3.8) 



where /3 <G [0,1]. 

Therefore, one ends up with the following scaling be- 
haviour for x , 



X 1 




+ VT- 7 ] 
/3)T~ 5 + VfiT] 
+ (1 - /3)T- 5 ] 



regime 1, 
regime 2, 
regime 3, 
regime 4, 
(3.9) 

where n p h — N p h/ L 2 and n k = N k /L are the phonon and 
kink densities respectively. The script letters indicate 
other quantities than the ones expressed already in the 
equations above. 



IV. COMPARISON WITH EXPERIMENTAL 
DATA AND CONCLUSION 

We compare our qualitative theoretical results with ex- 
perimental data in Ref. Q]. In figure 1 of this reference 
one sees that for a sample of highly purified lead which 
has been plastically stretched at low temperatures, the 
thermal conductivity at low temperatures has a pecu- 
liar shape: up to certain temperature it increases with 
temperature, then starts decreasing and for even higher 
temperatures it starts increasing with temperature again. 
Annealing can make this effect less pronounced, but it 
seems not to be able to completely remove this feature. 
Assuming that (3 is neither or 1 and taking numerical 
results into account^, one sees from Eq. (|3.9p that for low 
temperatures x scales as 



T 



n k C + T 4 ' 
for higher temperatures as 

n k C + T 4 : 
for even higher temperatures as 

T -i 



n k C + T- 2 ' 



(4.1) 



(4.2) 



(4.3) 



rp5 



n k C + T 4 ' 



(4.4) 



So at the highest and semilowest temperatures, the scal- 
ing behaviour is the same. The exact prefactors are dif- 
ferent of course. This mimics the behaviour shown in 
the experimental data. In the semi-highest temperature 
regime the thermal conductivity will decrease with tem- 
perature, while in the other regimes the thermal conduc- 
tivity will increase with temperature. 

When comparing curves 6 and 7 in figure one, one sees 
that curve 6 and 7 have similar behaviour for higher tem- 
peratures. For lower temperatures though, curve 6 lies 
under curve 7. As curve 6 shows the thermal conductiv- 
ity for a sample which has been deformed, while curve 7 
shows the thermal conductivity for a lead sample which 
has not been deformed at all, this is in full agreement 
with the theory. The power-law for the thermal conduc- 
tivity for a sample with none or very little kinks has a 
lower power than that for a sample with many kinks. 
Therefore it makes sense that for low temperature, the 
thermal conductivity for a sample with many kinks is 
lower than that for a sample with very little kinks. For 
this observation, we can therefore conclude that samples 
which have not been plastically deformed at all show a 
much weaker version of this effect, proving that this ef- 
fect is indeed caused by phonon-kink scattering. This 
also shows that only a small amount of kinks are needed 
to let this effect appear. 

The experimental data for the normal state does not 
match with our theoretical calculations at all, since in 
the normal state the phonon contribution to the heat 
flux transport is much weaker than the electron contri- 
bution. Therefore the effect of phonon-kink scattering is 
not visible in that case. 

We thus see that the results of our model qualitatively 
agree with the experimental data. For a quantitative 
comparison we refer to Ref. 0- The work of S.I. Mukhin 
is in part supported by RFFI grant 12-02-01018. The 
work of J. A.M. van Ostaay was supported by an ERC 
Advanced Investigator Grant. 
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